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Abstract. The aim of the paper is to derive the numerical least-squares esti- 
mator for mean and variance of random variable. In order to do so the following 
questions have to be answered: (i) what is the statistical model for the estima- 
tion procedure? (ii) what are the properties of the estimator, like optimality 
(in which class) or asymptotic properties? (iii) how does the estimator work 
in practice, how compared to competing estimators? 



1. Introduction 

The kriging model based on the stationary random process V = {Vj] Ni 9 j D i = 
1, . . . , n} with the estimation statistics Vj — UjVi with an unknown constant mean 
to and variance a 2 and some correlation function p for the asymptotic solution 

lim E{[(Vj - to) - (Vj - to)] 2 } = a 2 

has the well-known (co-ordinate independent) least-squares disjunction. The aim 
of the paper is to find (on computer) a co-ordinate dependent disjunction of kriging 
model for non-asymptotic solution 

E{[(V j -m)-(V j -m)] 2 }=a 2 . 
2. Correlation function estimator 



Following 



j(h) = -E{(Vj - V j+h y} = -E{Vf-2V j V j+h + Vf +h } 

= E{V 2 } — E{VjVj + h} 

= E{V 2 } - E 2 {V 3 ) - {E{V 3 V. J+h } - E 2 {V,}) 
= a 2 - C(h) > 

we get 

< C(h)/C(0) = C(h)/a 2 = 1 - j(h)/a 2 = \p(h)\ . 
Since the correlation function estimator must be non-increasing 

\p n (h)\ = l-%(h)/a 2 
then only non-decreasing outcomes of the experimental semi-variogram 

^ ^ n—h 

= 2 (n-h) ^ ^ ~ Vj+h ^ k = °' ' ' ' ' " ~ 1 

1 ' 3 = 1 
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should be taking into consideration for h < d 

Sr 2 = %(d) > %(d - 1) > %(d - 2) > . . . > 7„(1) > 7„(0) = 

Since 

C(h) = E{VjV i+h } - E{Vj}E{V j+h } 

then 

n — h ^ n — h n—h 

dn{h) = — ^< VjVj+h J^w £ Vj S w ' +fc 

and the correlation function estimator is also bounded by 

|p„(/0| = A\(o) fc = o,...,d. 



3. Central Limit Theorem 
Let us consider a set of frozen in time t stationary random processes 

V t = {Vj; Ni 9 j D i = 1, . . . , * e N n+ i 
since from the Central Limit Theorem holds 

[iY,v ] t=n+1 - E{V 3 } t=n+1 



U t =n 



+1 



y/E{V?} t=n+1 - ^{^} t =„+l 



JV(0;1) 



7/*= 



^E{V}} t = n+k - ^{^}t=n+fc 

an independent set of theirs outcomes 



JV(0;1) 



for 



t = n + 1, . . . ,n + fc 
also must follow the Standard Normal distribution. 

Example. Let us consider long-lived asymmetric index profile recorded by 359 
close quotes of Warsaw's Stock Market Index shown in Fig. 1 . Let the values Vj — 
w ii5+ij ■ • ■ 7 ^359 on the latter asymmetric side be supposed to be response averages 
[vj]t of frozen in time t stationary random processes based on Vi = v\, . . . , vus 

V t = {V j; Ni3jDt=l 115} t ; * G Nut+i 

with the correlation function 

' -l-i" 1 - 0135 ^/*] 2 , forA„>0, 



Pti^ij) = 



+1, 



for Aij = 0, 
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Figure 1. Long-lived asymmetric index profile, Warsaw's Stock 
Market Index (WIG) from 11 X 1999 up to 19 III 2001 (359 close 
quotes) the numerical least-squares estimator of mean at j > n + 1 
(black dots) compared for frozen model to the asymptotic gen- 
eralized least-squares estimator of mean (grey line). The dashed 
vertical line represents j = n = 115. 

with unknown mean m t and variance of solving on computer the least-squares 
constraint in one unknown j for every t 



E{[{Vi-m)-(Vi 
substituting the solution into 

Vj - ufa 



m)] 2 } = o 2 



\ju)v 2 - {w)Vi) 2 



'115 



we get, for every t = 115+1,... ,115 + 102, the Significance Level for computer 
K-S test 

a < 0.19 . 
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